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ABSTRACT 

Context. Studies of the properties of the inner Galactic Bulge depend strongly on the assumptions about the interstellar extinction. 
Most of the extinction maps available in the literature lack the information about the distance. 

Aims. We combine the observations with the Besancon model of the Galaxy to investigate the variations of extinction along different 
lines of sight towards the inner Galactic bulge as a function of distance. In addition we study the variations in the extinction law in 
the Bulge. 

Methods. We construct color-magnitude diagrams with the following sets of colors: H-Ks and J-Ks from the VVV catalogue as 
well as Ks-[3.6], Ks-[4.5], Ks-[5.8] and Ks-[8.0] from GLIMPSE-II catalogue matched with 2MASS. Using the newly derived 
temperature-color relation for M giants that match better the observed color-magnitude diagrams we then use the distance-color 
relations to derive the extinction as a function of distance. The observed colors are shifted to match the intrinsic colors in the Besancon 
model, as a function of distance, iteratively thereby creating an extinction map with three dimensions: two spatial and one distance 
dimension along each line of sight towards the bulge. 

Results. Colour excess maps are presented at a resolution of 15'x 15'for 6 different combinations of colors, in distance bins of 
1 kpc. The high resolution and depth of the photometry allows to derive extinction maps to 10 kpc distance and up to 35 magnitudes 
of extinction in Ay (3.5 mag in Ak s ). Integrated maps show the same dust features and consistent values with the other 2D maps. 
Starting from the color excess between the observations and the model we investigate the extinction law in near-infrared and its 
variation along different lines of sight. 

Key words. Galaxy: bulge, structure, stellar content - ISM: dust, extinction 



1. Introduction 

Interstellar extinction is a serious obstacle for the interpretation 
of stellar populations in the Galactic bulge. It shows a non ho- 
mogeneous clumpy distribution and is extremely high towards 
the Galactic Center. 

Several studies of cumulative extinction at the distance of the 
inner Galactic Bulge have been made: Schultheis et al. (1999) 
used the DENIS near-infrared data set in combination with the- 
oretical RGB/AGB isochrones from Bertelli et al. (1994). They 
showed that for an appropriate sampling area the observed se- 
quence matched well with the isochrone with appropriate red- 
dening. The maximum extinction that could be reliably derived 
from the J and Ks data available from DENIS was about A v = 
25 m . They found that in some areas, presumably with the high- 
est extinctions, there were no J band counterparts for sources 
detected in Ks. For these regions, only lower limits to the ex- 
tinctions could be obtained. Dutra et al. (2003) obtained similar 
results using the same technique with 2MASS data. Gonzalez 
et al. (2011) mapped the extinction along the minor axis of the 
Bulge based on the Vista Variables in the Via Lactea (VVV) ESO 
public survey data (Saito et al. 2012) using the red clump gi- 



ants. Gonzalez et al. (2012) presents the full extinction map over 
the entire VVV data set covering 315 sq. degrees using the red 
clump technique. This map extends beyond the inner Milky Way 
bulge, to also higher Galactic latitudes, and is the most complete 
study extending up to A v ~ 35 m , much higher extinction values 
than most previous studies. Recently, Nidever et al. (2012) deter- 
mined high-resolution A Ks maps using GLIMPSE-I, GLIMPSE- 
II and GLIMPSE-3d data based on the Rayleigh-Jeans Color 
Excess ("RJCE") method. 

However, most of the studies are restricted to 2D extinc- 
tion maps. Drimmel et al. (2003) built a theoretical large scale 
three dimensional Galactic dust extinction map. Marshall et al. 
(2006) presented a 3D model of the extinction properties of the 
Galaxy by using the 2MASS data and the stellar population syn- 
thesis model of the Galaxy, the so-called Besancon model of 
the Galaxy (Robin et al. 2003). However, their study is limited 
by the confusion limit of 2MASS in the galactic Bulge region 
and the limiting sensitivity of 2MASS in high extincted regions 
(A v > 30 m ). Recently, Robin et al. (2012) did a major improve- 
ment in the Besancon model by adding a bar component which 
results in an excellent agreement by comparing 2MASS color- 
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magnitude diagrams as well as in the metallicity and radial ve- 
locity distribution (see Babusiaux et al. 2010, Gonzalez et al. 
2012, Uttenthaler et al. 2012 ). We will use this updated model 
to build up a new 3D extinction model using the GLIMPSE- 
II data which maps the Galactic Plane within ±10 degrees in 
longitude from the Galactic Center with a wavelength coverage 
between 3.6 to 8.0 fj.m allowing to map the highest extincted re- 
gions as shown in Schultheis et al. (2009). In addition, we use the 
recently published VVV data with a pixel size nearly 10 times 
smaller and a sensitivity in Ks of 3-4 mag deeper than 2MASS 
to map the 3D extinction in the near IR. 

The infrared color-color diagram J-Ks vs. K-[mid-IR] has 
been used by Jiang et al. (2006) and Indebetouw et al. (2005) 
to determine the extinction coefficients A m id-i r along different 
lines of sight. Gao et al. (2009) traced the extinction coefficients 
based on the data from GLIMPSE survey and found systematic 
variations of the extinction coefficients with Galactic longitude 
which appears to correlate with the location of the spiral arms. 
Zasowski et al. (2009) found a Galactic radial dependence of 
the extinction law in the mid-IR using GLIMPSE data while 
Nishiyama et al. (2009) determined the extinction law close to 
the Galactic Center between 1.2 to 8.0 //m. 

In this paper, we first introduce our data set in Sect. 2. 
We discuss in Sect. 3 the Besancon stellar population synthe- 
sis model together with the recent improvements. In Sect. 4 we 
describe the method to determine the extinction coefficients and 
in Sect. 5 the method to derive the 3D-extinction. We finish with 
Sect. 6 with the determination of the extinction coefficient part 
as well as the conclusion (Sect. 7). 

2. The data set 

2.1. The GLIMPSE-II survey 

The Galactic Legacy Infrared Mid-Plane Survey Experiment 
(GLIMPSE-II) surveyed about 20 square degrees of the cen- 
tral region of the Galactic Inner Bulge using the Spitzer Space 
Telescope (Werner et al. 2004) equipped with the Infrared Array 
Camera (IRAQ (Fazio et al. 2004). It surveyed approximately 
220 square degrees of the Galactic plane at four IRAC bands 
[3.6], [4.5], [5.8], and [8.0], centered at approximately 3.6, 
4.5, 5.8 and 8.0 //m respectively. The GLIMPSE-II data, cov- 
ers a range of longitudes from -10° to 10° and a range of lat- 
itudes from +1 deg to +2 deg depending on the longitude. The 
GLIMPSE-II coverage does not include the Galactic center re- 
gion | / | < 1° and | b \ < 0.75°, which was observed by the 
GALCEN program (Ramirez et al. 2008). 

We use here the point source archive (GLMIIA, Churchwell 
et al. 2009). The GLMIIA catalogs consist of point sources with 
a signal to noise higher than 5 in at least one band and less strin- 
gent selection criteria than the Catalog (see Churchwell et al. 
2009 for more information). It contains the entire GLIMPSE 
II survey region including the GALCEN data and GLIMPSE 
I data at the boundary of the surveys (at 1=10° and 1=-10°). 
Moreover, these catalogs were also merged with 2MASS result- 
ing in a seven bands catalog: three 2MASS bands and four IRAC 
bands. The photometric uncertainty of the GLIMPSE-II data is 
typically better than 0.2 mag and the astrometric accuracy better 
than 0.3". 

2.2. The VVV survey data 

The other data set comes from the Vista Variables in the Via 
Lactea (VVV) ESO VISTA near-IR public survey. Observations 
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Fig. 1. Star counts of GLIMPSE II data of two example fields for the 
2MASS bands and IRAC bands. Each field is 15' x 15'. 
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Fig. 2. Star counts of VVV survey data of two example fields for the 
2MASS bands. Each field is 15' x 15'. 



were carried out using the VIRCAM camera on VISTA 4.1m 
telescope (Emerson & Sutherland 2010) located at ESO Cerro 
Paranal Observatory in Chile. The five year variability campaign 
of the survey, which started observations during 2010, will re- 
peatedly observe in Ks band an area of ~ 564 deg 2 , including a 
region of the Galactic plane from 295° < 1 < 350° to -2° < b < 
2°, and the Bulge spanning from -10° < 1 < 10° and -10° < b < 
+5°. During the first year of observations, the complete survey 
area was fully imaged in five bands: Z, Y, H, J, and Ks. In this 
work, we make use of the J, H and Ks catalogs for the same field 
of the GLIMPSE-II data from the 314 sq.deg coverage of the 
Galactic bulge. For a complete description of the survey and the 
observation strategy we refer to Minniti et al. (2010). A detailed 
description of the data products used in this work, corresponding 
to the first VVV data release, can be found in Saito et al. (2012). 

A description of the single-band photometric catalogs, pro- 
duced by the Cambridge Astronomical Survey Unit (CASU), 
and the procedure to build the multi-band catalogs can be found 
in Gonzalez et al. (201 1). These catalogs were used in Gonzalez 
et al. (2012) to construct the complete 2D extinction map of the 
Bulge and are the same ones used for the present study. We 
provide here only a brief summary of the procedure. J, H, and 
Ks catalogs produced by the CASU are first matched using the 
SLITS (Taylor 2006) code based on the proximity of sky posi- 
tions. Since VVV photometry is saturated for stars brighter than 
Ks=12 mag, 2MASS photometry is used to complete the bright 
end of the catalogs. Therefore, VVV catalogs are calibrated into 
the 2MASS photometric system. As described in Gonzalez et 
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al. (2011), this calibration zero point is obtained by matching 
VVV sources with those of 2MASS in the magnitude range of 
12 m < Ks < 13™. This magnitude range was selected to avoid 
the saturation limit of VVV and the faint 2MASS magnitudes 
with poor photometric quality. This last point is also reinforced 
by the usage of only sources flagged with high photometric qual- 
ity in the 2MASS catalogs. The calibration zero point is calcu- 
lated and applied individually for each VVV tile catalog. The 
final J, H, Ks magnitudes in the calibrated catalogs are in the 
usual 2MASS system. 

2.3. Completeness limit 

Due to the different pixel size of the different surveys 
(GLIMPSE-II, 2MASS, VVV), completeness limits changes as 
a function of wavelength and depends on the field direction. 
We divided the GLIMPSE II fields to subfields of 15' x 15' 
each to ensure a sufficient number of stars (in the order of at 
least 200) for our fitting process (see Sect. 4). The complete- 
ness limit for every sub-field was then calculated in each fil- 
ter. In order to estimate the completeness we use star counts as 
a function of magnitude where the peak of each histogram in 
each filter gives the corresponding completeness limit. The cata- 
logue becomes severely incomplete beyond the point where the 
star counts reach maximum. We are aware that this is a simple 
treatment of the completeness limit. However, we see here also 
clearly how strongly the completeness limit changes along dif- 
ferent lines of sight as demonstrated (Saito et al. 2012) using 
artificial sources. 

Fig. 1 shows the star counts as a function of magnitude 
for GLIMPSE-II and 2MASS catalogues for two subfields, one 
located at (l,b)=(0,+l) which is our reference field (the so- 
called C32 field) due to its very homogeneous and low extinc- 
tion (Omont et al. 1999, 2003). The second field is located at 
(l,b)=(0,-l) selected here for its higher and more clumpy extinc- 
tion (Ojha et al. 2003). Due to crowding and extinction the com- 
pleteness limit can change within one or two magnitudes in the 
infrared bands. Clearly seen is also the difference in the com- 
pleteness limit at [3.6] and [4.5] compared to [5.8] and [8.0] for 
both fields indicating the higher sensitivity of [3.6] and [4.5]. 

Figure 2 shows the same fields this time with star counts 
histograms from the VVV data. Note the "artificial" raise at 
Ks = 12 m which is the upper brightness limit of VVV. We de- 
cide to use in the following analysis the VVV data only from 
Ks > 12 m . Figure 2 demonstrates a clear gain in sensitivity 
of VVV in comparison with 2MASS near-IR photometry. As 
shown by Saito et al. (2012), the VVV data is confusion limited. 
However, due to not fully matched spatial scales of the two cat- 
alogues and the fact that the overlap between 2MASS and VVV 
is in the magnitude range where the completeness between the 
two catalogues is very different, we do not investigate this issue 
further (this is beyond the scope of this paper), but it should be 
kept in mind when assessing the results. 

We divided the VVV data in the same sub-fields as in the 
GLIMPSE-II data set. 



3. The Besancon model 

3. 1 . Main description 

The Besancon Galaxy model is a stellar population synthesis 
model of the Galaxy based on the scenario of formation and evo- 
lution of the Milky Way (Robin et al. 2003). It simulates the stel- 
lar contents of the Galaxy by using the four distinct stellar popu- 



lations: the thin disc, the thick disc, the bulge and the spheroid. It 
also takes into account the dark halo and a diffuse component of 
the interstellar medium. For each population, a star formation- 
rate history and initial mass function are assumed which allows 
to generate stellar catalogs for any given direction, and returns 
for each simulated star its magnitude, color, and distance as well 
as kinematics and other stellar parameters. The model, in its 
present form, can simulate observations in many combinations 
of different photometric systems from the UV to the Infrared. 
The photometry is computed using stellar atmosphere models 
(Basel 3. 1) from Lejeune et al. (1997) and Lejeune et al. (1998). 
However, while the Basel 3.1 stellar library is reliable for tem- 
peratures higher than 4000 K, atmospheric models become much 
more complex due to the appearance of molecules such as TiO, 
VO, H2O for cooler stars. As we have a significant contribution 
of M giants in our study we decided to adapt a more realistic T e ff 
vs. color relation extending to cooler objects of about 2500 K. 
This new relation is described in Sect. 3.2. 

Recently, a new Bulge model has been proposed by Robin 
et al. (2012) as the sum of two ellipsoids: a main component 
that is the standard boxy bulge/bar, which dominates the counts 
up to latitudes of about 5°, and a second ellipsoid that has a 
thicker structure seen mainly at higher latitudes. The resulting 
2MASS color-magnitude diagrams as well as star counts of this 
two component model clearly indicate the large improvement 
(see Robin et al. 2012). Uttenthaler et al. (2012) show also an 
excellent agreement in the kinematics based on radial velocity 
measurements. 

For more detailed information about the Besancon model we 
refer to Robin et al. (2003) and its most recent update in Robin 
et al. (2012). 
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Fig. 3. Comparison of the the simulated color-magnitude diagrams us- 
ing the old Tj. ff vs. color relation (left) compared to the new one (right) 
for the reference field (1=0.00 and b=1.00). The colorbar indicates the 
number of the stars. 
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Fig. 4. (a): log T c tf vs log g grid in the Besancon model. Black triangles 
indicate the coverage of the Besancon model (Robin et al. 2003) while 
red stars show the new extended grid, (b): J-Ks vs. T e jf relation. Black 
triangles show the old relation while red crosses the new one. The blue 
symbols is the relation for M giants from Houdashelt et al. (2000). (c). 
Similar as (b) but for K-[4.5] vs. T cff . 



3.2. Implementation of IRAC colors and a new 
temperature-color relation 

Reliable temperature-color relations are of extreme importance 
for comparison of the synthetic color-magnitude and color-color 
diagrams with the observed ones. We used the reference field 
at (l,b)=(0,+l) to check the new synthetic colors in the [IRAC] 
filters. Fig. 3 shows a kink in the Ks - [4.5] vs. [4.5] color mag- 
nitude diagram at [Ks] - [4.5] ~ +0.1, [4.5] ~ 10 which we 
identify clearly also in the corresponding temperature vs. color 
relation (see Fig. 4c). While the BaseO.l library predicts for 
T e ff < 3000 K decreasing J-Ks values, our new adopted re- 
lation shows a steady increase in J-Ks until T e s = 2500 K. 
We also noticed that the full parameter space in the logT e ff vs. 
log g plane was not covered in the previous model, especially 
there is a gap in the temperature range between 2500-4000 K 
for 0.5 <logg<3. 

We updated the temperature-color relation and filled in the 
gaps, by using the isochrones from Girardi et al. (2010). These 
isochrones include the TP-AGB phase (Marigo et al. 2008) in- 
cluding mass-loss. Starting from a group of isochrones with 
ages ranging between 5-10 Gyr we derived a new tempera- 
ture vs. color relation for metallicities between [Fe/H]=-2 and 
[Fe/H]=0.0. These ages and metallicities are characteristic of 
the Bulge populations (e.g. Zoccali et al. 2003, Zoccali et al. 
2008, Bensby et al. 2011). The comparison between the color- 
temperature relations from the Besancon model (so called "old" 
relation) and the one derived by us (so called "new" relation) is 
shown in Fig. 4b and c. For comparison we show the J-Ks vs. 



T e ff relation obtained by Houdashelt et al. (2000) for M giants. 
The clear difference between the new and old temperature vs. 
color relations for T e ff < 4000 K is obviously a culprit for the 
bent in the color-magnitude diagram (Fig. 3), which disappears 
after adopting our new relation. We also see that using the new 
relation the J-Ks colors will be systematically redder by about 
0.05 mag. 

Figure 5 shows the comparison of a color-magnitude dia- 
gram in a low-extinction field -the so-called Baade's window. 
We use here the catalog of Uttenthaler et al. (2010). One notices 
clearly that the predicted color in the "old" model is too blue 
compared to the observations and we see in addition a much 
larger fraction of foreground objects (Ks - [4.5] < 0.5) which 
we don't find in the data. 

4. Method 

We obtain the extinction by using the color excess between the 
intrinsic color (derived from the model defined as Coj ns ) and the 
observational apparent color (derived from the observation, here 
we defined as Co t, s ). Assuming an extinction law, we get: 

A A = c A x (Co obs - Co ins ) (1) 

where q is related to the extinction coefficient. As we trace 
our fields close to the Galactic Center, we used the extinction 
coefficients from Nishiyama et al. (2009). We ignore here the 
effect of the broadband extinctions and the reddenings on the 
SED of the individual sources as our method does not allow for 
a more sophisticated treatment which could imply systematic 
errors. This results in: 

A Ks ,j- Ks = 0.528 x ((J - Ks) ohs -(J- Ks) ins ) 
A Ks ,h-ks = 1.61 X ((H - Ks) obs -(H- Ks) ins ) 
A [3 .6] M .6] = 1-005 x ((Ks - [3.6]) obs - (Ks - [3.6]) ins ) 
A [4 .5],^-[4.5] = 0.640 x ((Ks - [4.5]) obs - (Ks - [4.5])^) 
A [5 . 8 ], M 8] = 0.562 x ((Ks - [5.8]) o6s - (Ks - [5.8])^) 
A m , Ks - [m = 0.748 x ((Ks - [S.0])obs ~ (Ks - [8.0])*,) 

In the following we apply the same method as Marshall 
et al. (2006) assuming that the distance grows as the appar- 
ent color/extinction increases. Fig. 6b,c shows an example of 
this relation for our reference field. Contrary to Marshall et al. 
(2006) we use the full information of the stellar populations in 
the CMD, e.g. we include also the M dwarf population. While 
the local dwarf population (d < 1-2 kpc) can be easily identi- 
fied (see Marshall et al. 2006) by simple photometric color cri- 
teria (e.g. J-Ks), this is not the case for dwarfs located at larger 
distances where they cannot be separated from giants. This is 
especially true for the VVV data where we are able to detect 
dwarf stars up to larger distances. However, the model simula- 
tions show that the contribution of dwarf stars within the VVV 
completeness limit used in our study is rather small (~ 1 %) but 
can increase up to 20-30 % going to fainter magnitudes beyond 
the VVV completeness limit. 

The method can be described in following steps: 

1. First we add photometric errors and a diffuse extinction to 
the intrinsic colors Co,„ s from the model for each subfield. 
The assumption of a diffuse extinction is necessary to ensure 
that the colors of the model increase with distance. We as- 
sume a diffuse extinction of 0.7 mag/kpc in the V band. We 
use exponential errors for the GLIMPSE-II data in all seven 
filters. Fig. 6a shows an example of the photometric errors 
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Table 1. Extinction as a function of Galactic longitude and latitude as well as distance based on VVV data. 

1 b E(J - Ks)j_i 0kpc E(H — Ks)i-\Q klK crE(J - Ks)i- iakpc crE(H - Ks)i_ lokpt . 
Notes. For each position we give the E(J-Ks), E(H-Ks) as well as the corresponding sigma for each distance bin starting from 1 to lOkpc. 

Table 2. Extinction as a function of Galactic longitude and latitude as well as distance based on GLIMPSE-II data. 

1 b £(jtj-[3.6])i_i 0tpc E(Ks-[4.5]),_ IQt „ c E(Ks - [5.8])i_i 01) ». E(K S - [8.0])i_i, u ,,. uE(K - [3.6])i-nu )x crE(K S - [4.5])^ lMpc aE(Ks - Pggjpjg <jE(Ks - [8.0])i_i 0t)x . 

Notes. For each position we give the E(Ks-[3.6]), E(Ks-[4.5]), E(Ks-[5.8]), E(Ks-[8.0]) as well as the corresponding sigma for each distance bin 
starting from 1 to lOkpc. 
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as a function of the [3.6] magnitude for our reference field 
(l,b) = (0.00, +1.00). For the VVV data, the error is much 
smaller and we adopted a constant error of 0.05 mag (Saito 
et al. 2012). This procedure provides a first set of simulated 
colors (Co simiini ) 

Colors are then sorted, both for observed and simulated 
data (Co bs and Co s i m i n i), and normalized by the corre- 
sponding number of colors in each bin given by «/,,„ = 
floor(min([N o b s ,Nbes])/100). N b s and Nbes are the total 
number of observed stars (in each subfield) and of stars in 
the Besangon model, respectively. This ensures that we have 
at least 100 stars per bin. 

In each corresponding color bin, the Co lns of the stars from 
the model data and the Co„b s from the observational data are 
used to calculate the extinction using Eq. 1, as well as the 
median distance from the model. A first distance and extinc- 
tion estimate can be thus obtained (A Is ,). Due to the satura- 
tion limit of VVV (Ks = 12), we do not get enough stars in 
the first distance bins (0-3 kpc). Therefore our 3D map starts 
at 4 kpc for VVV, while at 3 kpc for GLIMPSE-II data. 
The extinction as well as the distance estimated from the 
third step is then directly applied to the intrinsic magni- 
tudes in the model which gives us a new simulated color 
{Co S i m> \ = Co ins + A\ st ). 

We construct histograms for each color of the new simulated 
data (Co sjm {) and the observational data (Co ins ) using the 
same binsize (0.05 mag). The^f 2 statistics (Press et al. 1992; 
Marshall et al. 2006) are then used to evaluate the similarity 
of both histograms, given by: 



X 1 = Yj { ^ N ° bs!Ns > 



's/Ns im IN obsH-obsj) /(.^sim: + Mobs 



where N b s and N S i m are the total number of observed and 
simulated stars in each subfield, while n b Sj and n sim . show 
the number of stars for the jth color bin of the observations 
and simulated data, respectively. 

With the first extinction estimate applied we have now a new 
set of simulated data. For this new set of data, some of the 
stars may be fainter to fall outside the completeness limit 
while others of them may be brighter and come inside the 
limit. We go back to the second step with the new Co s ; m ,i 
which would be used iteratively to refine our results. In total 
we perform 20 iterations for each sub-field and each color. 
The minimum x 2 is men our final result. We have chosen 
20 loops to ensure that we always find a minimum value. In 
most of the cases the x 1 converges after one or two loops to 
a minimum and then becomes stable in the following loops. 

Fig. 6c shows an example of the iteration process using the 
Ks-[3.6] color of the GLIMPSE-II data in our reference field 
(l,b)=(0.00, 1.00). We see clearly the difference between the first 
extinction estimate and our final extinction determination after 
20 loops. Note that the we adapt the same completeness limit as 
in the observations for the model simulations (see Sect. 2. 3). 

To test our method, we used artificial data by adding manu- 
ally extinction values to the intrinsic magnitude of the model. 
We used different lines of sight with a different number of 
clouds. Figure 7 shows two examples of several clouds. For 
field (l,b) at (0.00, 1.00), one located at a distance of 2 kpc with 
Ay = 4.5 mag, the second one at 6 kpc with Ay = 2.5 mag. For 
the field (l,b)=(5.00,0.00), we assumed three clouds: A v = 2.20 
mag at 2.0 kpc, 2.50 mag at 5.0 kpc and 2.70 mag at 7.0 kpc. 
)(;We applied our method to our simulated data. We got 328.000 
sources for the first field and 182.000 for the second field. We 
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Fig. 6. The process to calculate the extinction of the example subfield 
(l,b) of (0.0,+ 1.0). From upper left to lower right, they are:(a) the er- 
rors of [3.6] vs. the [3.6] diagram, (b) the color Ks-[3.6] vs. distance 
diagram, (c) the extinction at [3.6] vs. the distance diagram where we 
show the first initial guess of extinction after the first iteration (black) 
and the final iteration after 20 loops (red) and (d) the color Ks-[3.6] 
distribution of the data compared to the first guess iteration (black) and 
the final one (red). 
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Fig. 7. The test of our method. The red line is the artificial extinction 
with clouds. The black crosses show our derived extinction values. The 
left panel is of field (l,b)=(0.00,+1.00) and the right (l,b)=(5.00,0.00). 



cut the data at Ks = 18 mag artificially, which goes much deeper 
than the observational data. Figure 7 shows that our method suc- 
cessfully derives the extinction along the distance and that we 
are able to trace this multi-cloud structure. Modifying the input 
parameters of the galactic model such as galactic scale length, 
the size of the hole of the galactic disc or changing the Bulge lu- 
minosity function could lead to systematic differences in our de- 
rived extinction map. But this systematics, as shown in Marshall 
et al. (2006), is of the order of the random uncertainty of the 
method. 



5. Results and discussion 

The full results are listed in Table 1 and Table 2, both avail- 
able in electronic forms at the CDS 1 . Each row of Table 1 and 
Table 2 contains the information for one line of sight: Galactic 
coordinates along with the measured quantities for each distance 
bin E(J-Ks), E(H-Ks) and respective uncertainties (VVV, Table 
1) as well as E(Ks-[3.6]), E(Ks-[4.5]), E(Ks-[5.8]), E(Ks- 
[8.0]) and their uncertainties (GLIMPSE-II, Table 2). These re- 
sults will be also added into the BEAM calculator 2 webpage 
(Gonzalez et al. 2012) where the 2D VVV extinction map is al- 
ready available to the community. Users of the BEAM calculator 
can choose to retrieve the extinction calculation adopting a spe- 
cific reddening law and distance. 

5.1. Extinction along different lines of sight 

As an example of our results we discuss in more detail here two 
specific lines of sight: one the reference field (l,b)=(0.00,1.00) 
and the other is the field at (l,b)=(0.00,-1.00). 

Figure 8 and Figure 9 show the corresponding color- 
magnitude diagrams and color distribution for the observational 
data and the simulated data with the best-fitted final extinction 
in each filter. Besides the general good agreement in the shape 
of the CMD and in the color distribution, we want to emphasize 
also that the number densities in the CMDs between the obser- 
vations and the model agree very well. However we notice some 
interesting features by comparing the CMDs between the model 
and the observations: 

- We see in the [5.8] vs. Ks-[5.8] and especially in the [8.0] vs. 
Ks-[8.0] GLIMPSE-II diagram clearly the Asymptotic Giant 
Branch phase. The tip of the AGB is approximately at [8.0] 
a bit brighter than 8.0 mag. This phase is not reproduced by 
the model and needs to be improved. 

- Contrary to GLIMPSE-II, the model predicts too many stars 
for VVV in all filters at the fainter end of the CMD. We sus- 
pect that the completeness limit of VVV is underestimated. 
Saito et al. (2012) did for two VVV fields a complete analy- 
sis of the completeness limit by using artificial sources. They 
found that the source detection efficiency reaches 50% for 
stars with 16.4 < Ks < 16.9 mag for the tile b314 located at 
(l,b)=(352.21,-0.92). Our derived completeness limit is Ks = 
14.9 which reaches according to Saito et al. (2012) the 75% 
completeness level. A detailed analysis of the completeness 
limits of the VVV survey using artificial sources would be 
needed. 

- We notice that the model has produced a second branch in 
the J-Ks vs Ks diagram at about Ks = 12.5 mag which we 
do not find in the VVV data (see Fig. 9). This branch is re- 
lated mainly to the galactic disc K-giant population, which is 
clearly too prominent in the model. In contrary to the J-Ks 
and H-Ks color, the J-H color shows a large scatter in the 
temperature vs. color relation. This scatter is mainly due to 
the large influence of J-H to metallicity and gravity (see also 
Bessell et al. 1989). We decided therefore not to use the J-H 
color in our analysis. 

Figure 10 to Figure 13 show how interstellar extinction 
changes with distance in the different bands for four subfields. 



1 Table 1 and Table 2 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/. 

2 http://mill.astro.puc.cl/BEAM/calculator.php 
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Fig. 8. CMD and color distribution for all the colors of the field 1=0.00, b=1.00. The lower panel shows the observed CMD, the middle panel 
the synthetic CMD with the best-fitted extinction and the upper panel the histograms in the color distribution. The red color scale (lower pannel) 
shows the observations (noted as ' VVV or GLIMPSE' in the title of the inner colorbar) while the blue for simulated data (noted as 'model'). The 
dashed lines show the completeness limit. The colorbar inside the CMD shows the number of stars of each corresponding dataset. 
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Fig. 9. CMD and color distribution for all the colors of the field 1=0.00, b=-1.00, same as Fig. 8. 
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Fig. 11. The distance vs. extinction diagram for sub-field 1=0.00, b=- 
1.00. The symbols are the same as in Fig. 10. 



Fig. 13. The distance vs. extinction diagram for sub-field 1=0.00, 
b=+1.75. The symbols are the same as in Fig. 10. 



(Wall & Jenkins 2003). We re-sampled the original datasets of 
the observational and the model data 1000 times. The r.m.s scat- 
ter gives the uncertainty of our results. Fig. 14 shows a typical 
example. It shows the distributions of our results obtained for 
the 1000 samples obtained for one colour bin (Step 3 in Sect. 4) 
for the GLIMPSE-II data in [3.6], left for the distance and right 
for the extinction. A Gaussian fit was made to the distributions, 
where the full-width-half maximum of the distribution gives the 
uncertainty in distance and extinction. This process is applied 
for each bin and each filter. The extinction changes along with 
the distance. As expected, it is higher for increasing distance. As 
seen in Fig. 10 to Figure 13, extinction gets flatter beyond 7 kpc. 
As shown by Marshall et al. (2006) this is probably due to the 
fact that inside the molecular ring the density of interstellar mat- 
ter is lower than outside. It has been shown also that this region 
contains the dust lanes of the bar (Marshall et al. 2009), which 
can explain that at positive longitudes one can see an increase of 
the absorbing matter after this plateau, while it appears further 
away at negative longitudes. 




6789 10 2 4 6 8 
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Fig. 12. The distance vs. extinction diagram for sub-field 1=10.00, 
b=1.00. The symbols are the same as in Fig. 10. 



Different symbols stand for different data set: squares for VVV 
data and triangles for GLIMPSE-II data while different colors 
denote the various filters. The errors of the derived distances 
and the extinction were calculated using the bootstrap method 




Fig. 14. The distribution of the distance (left) and the extinction (right) 
derived from 1000 samples with bootstrap method. 



As already mentioned in Sect. 2.3 the VVV data used here 
are only for Ks fainter than 12 mag. For this reason, the first bin 
of the VVV data is 4 kpc. The VVV data on the other hand goes 
much deeper (d ~ 10-12 kpc). But 2MASS starts to be incom- 
plete at 12 mag depending on the line of sight. Therefore, there 
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Fig. 15. 3d extinction and error maps for Ak s - X-axis denote galactic longitude while y-axis galactic latitude. The units are in <5Ak s kpc . The map 
shows the weighted average extinction converted from the seven individual colors excesses using the reddening law derived from Sect. 6. 



is in general a very small overlap between 2MASS and VVV 
which is also seen in the CMDs (see Figs. 8 and 9) where we 
miss in the VVV data the bright foreground branch of K giants 
clearly visible in GLIMPSE-II. For the four [IRAC] bands we 
see very similar extinction vs. distance relations. 

5.2. The three dimensional extinction map 

In order to visualize the distribution of extinction in 3D we di- 
vided the extinction between subsequent bins by the distance be- 
tween them similar as done by Marshall et al. (2006). The dis- 
tance intervals are interpolated in bins of 1 kpc. All the 3D-maps 



using the six different color excesses i.e. E(J - Ks), E(H - Ks), 
E(K- [3.6]), E(K- [4.5]), E(K- [5.8]), E(K- [8.0]) are presented 
in the Appendix A. The units are in kpc 1 . In addition we also 
show the error maps for each pair of filters. These errors are the 
combinations of the standard deviations of the observed colors 
and the synthetic ones. Note that on average the errors are lowest 
between 4-8 kpc. We see a very similar picture of the dust dis- 
tribution very concentrated to the galactic plane from all maps. 
As expected, extinction becomes higher as the distance from the 
Sun increases. 

We converted the seven color excess to Ak s using the ex- 
tinction law derived from Sect. 6 separately. The extinction map 
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Fig. 16. Dust extinction in the Galactic Plane for |b| < 0.25. The dotted lines correspond to distances of 5, 10 and 15 kpc. The Sun position is 
marked as "X" and the GC as "+". Indicated on the y-axis is the galactic longitude. The upper panel shows our results while the lower one those 
from Marshall et al. (2006). 
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Fig. 17. Comparison of derived color excess values from J-Ks (VVV) to 
H-Ks (VVV), Ks-[3.6] (GLIMPSE-II), Ks-[4.5] (GLIMPSE-II), Ks- 
[5.8] (GLIMPSE-II) and Ks-[8.0] (GLIMPSE-II). 

of the weighted mean values is shown in Fig. 15. The maps are 
presented as extinction in Ks band and are also in units of kpc" 1 . 
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Fig. 18. 3d extinction comparison with Marshall et al. (2006) and 
Drimmel et al. (2003), where the black solid line is for Drimmel et al. 
(2003) , black dashed line for Marshall et al. (2006) and the blue line 
is the result of our determination. The star symbols denotes the cal- 
culated distance of the corresponding red-clump position of Gonzalez 
etal. (2012). 



In order to visualize better the dust extinction along the line of 
sight, Fig. 16 shows a view of our results (upper panel) together 
with that from Marshall et al. (2006) (lower panel) from the 
North Galactic pole towards the galactic plane at |b| < 0.25 sim- 
ilar to Fig. 9 of Marshall et al. (2006). Our results show globally 
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similar structures compared to Marshall et al. (2006). Except for 
1 = 0° and 1 = -7° we notice beyond 6 kpc a very smooth and 
featureless dust extinction. Clearly seen in this comparison is 
also the much lower extinction until 6 kpc compared to Marshall 
et al. (2006) which will be discussed in Sect. 5.3. The highest 
concentration of dust is around the Central Molecular zone, a gi- 
ant molecular cloud complex with an asymmetric distribution of 
molecular gas (see e.g. Morris & Serabyn 1996). The dust ex- 
tinction map shows very similar non-axisymmetric structure as 
observed in OH and CO. However, the spatial resolution of our 
map is not high enough to trace features like the lOOpc ellipti- 
cal and twisted ring of cold dust discovered by Molinari et al. 
(2011). 

The highest extincted regions are best traced by E(H - Ks) 
using the VVV data as well as E(Ks - [3.6]) going up to Ar s = 
3.5 mag. Due to the high sensitivity of VVV, we are able to 
trace in 3D high extincted regions (such as star forming regions, 
molecular clouds,etc.) until distances of 10 kpc. 

Concerning the E(Ks - [8.0]) map we want to stress that at 
8 //m PAHs become an important component of emission which 
weakens the stellar flux (see e.g. Cotera et al. 2006). We get 
systematically higher extinction due to the PAHs emission. 

Fig. 17 shows the comparison in the two-dimensional space 
between the color excess using E(J - Ks) and the E(H - Ks), 
E(Ks - [3.6]), E(Ks - [4.5]), E(Ks - [5.8]) and E(Ks - [8.0]) to 
see if we introduce any systematic offsets between the different 
filter pairs. The points with bootstrapping errors larger than la- 
have been excluded. In general we see a very good relation with 
a difference not larger than 0. 1 mag. The dispersion is higher 
for the GLIMPSE-II data which is due to the larger photometric 
errors. 

5.3. Comparison with other 3D-maps 

We compare our results with that from Marshall et al. (2006) 
and Drimmel et al. (2003) for the extinction in the Ks band. 
Fig. 18 shows the extinction v.s. distance diagram for two dif- 
ferent fields, one located at (l,b)=(0.00,1.00) and the other at ex- 
tremely high extinction located at the Galactic Center (1=0.00, 
b=0.00) compared to the dust model of Drimmel et al. (2003) 
and to Marshall et al. (2006). 

We can see that they show the same trend as mentioned in 
Sect. 5.1. The black solid line is from Drimmel et al. (2003) and 
the dotted line from Marshall et al. (2006). In order to compare 
with Marshall et al. (2006), we converted their A Ks values as- 
suming the extinction coefficients of Nishiyama et al. (2009). 
Figure 18 shows that for our reference field (l,b)=(0,+l) we get 
a similar Ak s vs. distance relation for d > 6 kpc while Marshall 
et al. predicts a rather steep slope for d < 6 kpc. For the high ex- 
tinction field at (l,b)=(0,0) there is a quite significant difference. 
We obtain a steeper slope and we get systematically higher ex- 
tinction values for d > 4 kpc. We added also for comparison a 
data point of the 2D map of Gonzalez et al. (2012) on the di- 
agram marked as star symbol. The distance was calculated by 
taking the difference of the dereddened mean magnitude of the 
red clump and the absolute magnitude of the red clump which 
we assume to be Mk s = -1.55 mag. As one can see the esti- 
mated distances from Gonzalez et al. (2012) agree remarkably 
well with our distances. Schodel et al. (2010) used adaptive op- 
tics observations of the central parsec of the Galactic center in 
the Ks band to map the extinction. They obtained a mean extinc- 
tion value around Sagittarius A* of Ak s = 2.46 m with a distance 
to the GC of R = 8.02 kpc, while Fritz et al. (2011) derived 



A Ks = 2.42 m . Our derived value of A Ks = 2.34 m is within the 
errors of the cited papers. 

For low extinction fields, our results are in agreement with 
Marshall et al. (2006). We get only for larger distances slightly 
higher A Ks values than Marshall et al. (2006). However, for the 
high extincted fields such as the GC (right panel of Fig. 18), 
we notice rather significant differences: While the extinction 
seems to flatten at a distance of about 6 kpc for Marshall et al. 
(2006), we notice a steady increase in extinction reaching val- 
ues of A Ks = 2.5 m compared to only A Ks = 1.5 m for Marshall 
et al. (2006). This low limit is due to the large incompleteness 
of 2MASS in this high extincted region compared to VVV. We 
emphasize that compared to Marshall et al. (2006) our distance 
steps shown here are much smaller and we trace thus better the 
derived distance-color relations. 

Thanks to its deeper photometry and higher spatial resolu- 
tion, the VVV survey catalogue has stars at larger distances than 
2MASS. For a red clump star with M Ks = -1 .65 m as an example, 
a brightness at Ks = 12 m , the reliable measurement by 2MASS 
towards crowded regions, means a distance of about 3.5 kpc with 
AKs=lmag, and 5.5 kpc if no extinction is applied. Since the 
Marshall results are derived based on the 2MASS data, there are 
difficulties in probing far and highly extincted regions. While at 
Ks=15 mag in the case of the VVV survey, the red clump star 
can be detected at a distance of 13.4 kpc if Ak s =1 mag and 5.4 
kpc if Aks=3 mag. Considering the RGB and AGB stars that are 
mostly brighter than red clump stars, the VVV survey detected 
many stars at a distance further than 5 kpc. This depth of VVV 
makes the results at larger distance more solid. This also ex- 
plains the systematically higher extinction we obtained for the 
GC field at 1=0.0 and b=0.0 (see Fig. 18) compared to Marshall 
et al. (2006). 

Our results are clearly very different from those of Drimmel 
et al. (2003), in particular for the Galactic Center region. The 
Drimmel model predicts much higher extinction than ours as 
well as Marshall's. The dust model of Drimmel et al. (2003) 
depends very much on the dust temperature. However, in the 
Galactic Center region the gas pressure and temperature is higher 
than in the galactic disk (see Serabyn & Morris 1996) which ex- 
plains the large difference in the extinction values compared to 
Drimmel et al. (2003). Our map is not sensitive to dust tempera- 
ture like theirs but to the modeled K/M giant stellar population. 
As our map is restricted to a spatial resolution of 15'x 15', small- 
scale variations such as seen e.g. by Gosling et al. (2008) cannot 
be resolved. 



5.4. Two dimensional extinction maps 

We compare our extinction maps with the 2D-maps of Schultheis 
et al. (1999) and Gonzalez et al. (2012). While Schultheis et al. 
(1999) used mainly the RGB/AGB star population, Gonzalez 
et al. (2012) restricted themselves to the red clump stars to trace 
the extinction. For that purpose we integrate our extinction map 
up to a distance to 8 kpc. The resulting 2D-maps of the com- 
bined A Ks are shown in Fig. 19. The individual maps of the 
different color excesses are presented in the Appendix B. They 
show clearly very similar dust features as described in Schultheis 
et al. (1999) and Gonzalez et al. (2012). We see clearly the small 
scale variation of the interstellar dust clouds concentrated to- 
wards the galactic plane which is not symmetric. Besides the 
clear dust features of the Central Molecular Zone,we also see 
a clear correlation between the high extinction and the number 
of detected infrared dark clouds based on molecular line obser- 
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Fig. 19. 2d extinction and error map for Ak s integrated to a distance of 8 kpc. The map shows the weighted average extinction converted from the 
seven individual colors excesses using the reddening law derived from Sect. 6. X-axis denote galactic longitude while y-axis galactic latitude. 
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Fig. 20. Comparison of A Ks with Schultheis et al. ( 1 999) (top panel) and 
Gonzalez et al. (2012) (lower panel). 



vations (Jackson et al. 2008), especially at the region around 
(l,b)=(353.0, 0.0) with a large concentration of YSOs. 

For a more quantitative comparison we smoothed the 
Gonzalez et al. (2012) maps to the same spatial resolution as 
ours (15'x 15') and corrected for differences in the extinc- 
tion law. Figure 20 shows the comparison. We see an excel- 
lent agreement between these maps with no systematic offset. 
For Aks > 1 .5 mag the dispersion between our map and that of 
Schultheis et al. (1999) gets larger which is caused mainly by the 
limited sensitivity of the DENIS survey in the J and Ks bands. 
In contrast, the dispersion compared to the VVV extinction is 
rather small, even for higher Ak- We are thus confident that our 
extinction determination is highly reliable. 

Figure. 21 shows that the intrinsic dispersion grows with 
increasing extinction. Lada et al. (1994) demonstrated that the 
form of the observed cr^sp versus Ay relation can be used to 
place constraints on the nature of the spatial distribution of ex- 
tinction. However, photometric uncertainties as well as uncer- 
tainties in the extinction law dominate this relation. Note that 
the dispersion in VVV is much smaller than for 2MASS due to 
its much smaller photometric errors. A higher spatial resolution 
would be needed to study this feature more in detail. The 2D 
maps of Gonzalez et al. (2012) have 2'x 2'resolution resolution 
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Fig. 21. The errors v.s. the extinction diagram for VVV (left panel) and 
2MASS+GLIMPSE (right panel). The red solid line is an exponential 
fit. 
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but their dispersion is dominated by the distribution of dust along 
the distance in that line of sight. 

Figure. 22 shows the total number of observed stars for each 
V V V-subfield and the corresponding Ak s . We see clearly that the 
number of stars gives an indication of the extinction, i.e., lower 
numbers means higher extinction. 



6. The extinction coefficients 

As we determine the extinction maps using E(A-Ks) we can cal- 
culate the extinction coefficients. For the extinction coefficients 
A^/Afcs, we used here the values from Nishiyama et al. (2009) as 
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the initial values. If A is the wavelength we are going to derive 
and a is the assumed wavelength, we can write: 




10 

A [fj,m] 



Fig. 23. Our derived interstellar extinction curve together with the com- 
parison of different literature work. Note that we use the mean extinc- 
tion values over all fields. The red dashed line indicates a simple power 
law A A oc A~ [ 75 from Draine (1989) while the blue the power law of 
Cardelli et al. (1989) with A A oc / T L61 and the green the power law of 
Fitzpatrick & Massa (2009) with A A oc A~ iM . The dot-dashed line in 
magenta shows the extinction curve of Fritz et al. (201 1) in the GC. 



explain the flattening of the extinction curve between 3-8 /jm. 
As shown by Gao et al. (2009) (see their Fig. 7) a larger mean 
dust grain size (e.g. 0.3 //m) would give a steeper extinction law 
(i.e., smaller A^/Ak, ratios) than that of the diffuse ISM where 
the mean dust size is 0. 1 fim. Our derived extinction curve is very 
similar to that of Nishiyama et al. (2009). Interestingly, we find a 
smaller A[g.8]/AKs value than Nishiyama et al. (2009). However, 
our value is within the errors close to the value of Nishiyama 
et al. (2009). Fritz et al. (2011) studied the near-IR and mid-IR 
extinction law towards SgrA* using hydrogen emission lines. 
They found a power-law slope of alpha = -2.11 + 0.06 short- 
ward of 2.8 fim. This is in agreement with our derived value 
of a = -2.175 and those found in the literature. Figure 23 
shows that our derived extinction curve agrees remarkably well 
to that of Fritz et al. (2011) (indicated by the dot-dashed line) 
and confirms our low extinction value at [5.8]. As pointed out 
by Fritz et al. (2011) classical grain models (Li & Draine 2001, 
Weingartner & Draine 2001). using mainly silicate and graphite 
grains with different size distributions fail to reproduce the ob- 
served extinction law and especially the change of slope in the 
extinction law. Water ice features could be responsible for the 
slope change in the mid-IR. However further modeling is neces- 
sary. 

The comparison shows that our mean values in the IRAC 
bands are comparable to those found in the literature. However 
it depends on the assumed A a /A Ks . It can be seen from Table 3 
that A,i/Aks gets smaller value when A a /A^ s is bigger, which 
is a natural result from equation (3) because (A a / A^ s -l) is neg- 
ative for any wavelength longer than Ks. The values of Aj/A Ks 
assumed or derived have a minimum of 2.56, a maximum of 
2.89 and a mean of 2.70. In comparison, Gao et al. (2009) used 
Aj/A Ks =2.52. 



A A /A Ks = 1 + (AJA Ks - 1) x -i — i (3) 

E(a - Ks) 

A ff /A Ks is the initial value for the assumed wavelength from 
Nishiyama et al. (2009). E(/l-Ks) /(E(a-Ks)) is the color excess 
derived from the iteration of^- 2 test (see Sect.4). By Eq.(3), we 
make the assumed wavelength as J, H, [3.6], [4.5], [5. 8], [8.0] to 
calculate the extinction coefficients for each sub-fields. The val- 
ues we adopted are from Nishiyama et al. (2009) because their 
line of sight of study is towards the Galactic center. However, 
considering the variation of the extinction law which are sug- 
gested by a couple of studies, using a unique value of A a /Aa s 
could introduce some errors. 

6.1. The mean extinction coefficients 

Different choices of the assumed A a /A Ks would get different re- 
sults for Ax/Aks- We list our results in Table 3 varying the as- 
sumed extinction coefficient Aa/A^s- The mean extinction coef- 
ficients of all our fields are also listed in Table 3 and compared to 
Gao et al. (2009), Indebetouw et al. (2005) and Nishiyama et al. 
(2009). The assumed values are noted as bold fonts in Tab. 3. 
Figure 23 shows the derived extinction curve. The black circles 
indicate our mean value (see Tab. 3). The red line indicates a 
simple power law A A oc A~ l 15 from Draine (1989) while the blue 
line the power law of Cardelli et al. (1989) with A A oc A~ 1M 
and the green line the power law of Fitzpatrick & Massa (2009) 
with A A oc /T 1 84 . It is clear that these simple power laws can not 



6.2. Variation of the extinction coefficients 

The extinction coefficient may differ along different lines of 
sight. Gao et al. (2009) showed an interesting distribution of 
extinction ratios as a function of galactic longitude with 131 
GLIMPSE fields. In this work, we take the extinction coeffi- 
cients A a /Aks of all our sub-fields in our field. We take the 
median value of the A A /A Ks in each bin of longitude or lati- 
tude. To see the variation clearly, we present here the relative 
value (^). Figure 24b shows a 2D map of Ap.ej/A^ assum- 
ing Aj/A Ks = 2.86, and Fig. 24a,c show the variation of the 
extinction coefficients (i.e. A^^/Ak^A^^/Ak^A^^/Aks and 
A[8.o]/Aat. s ) as a function of galactic longitude (Fig. 24a) and lat- 
itude (Fig. 24c). Overplotted in figure 24a are also the values 
from Gao et al. (2009). Note that they use larger binsizes (1 de- 
gree) compared to ours. We notice in all [IRAC] bands a similar 
behavior in the variation of the extinction coefficient along the 
galactic longitude which follows the extinction curve variation 
of Gao et al. (2009). While this variation is rather small, we see 
a surprisingly large peak visible towards the galactic plane indi- 
cating a large variation with galactic latitude (see Fig. 24b). This 
variation seems to be higher for longer wavelengths except for 
8.0 /im where the amplitude becomes small again. This might 
be related to the 9.7 /mi silicate absorption feature (Gao et al. 
2009). Figure 24 indicates that the extinction curve becomes 
flatter in the mid-IR approaching the galactic Center region. A 
detailed comparison with interstellar dust models is needed to 
explain this feature which is beyond the scope of this paper. 
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Table 3. Extinction coefficients A A /A Ks . 



a 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


A([im) 


1.24* 


1.664* 


2.164* 


3.545 


4.442 


5.675 


7.760 


VVV 


fixed J 


2.86 ±0.01 


1.68 ±0.03 


1.00 










fixed H 


2.63 +0.18 


1.60 ±0.01 


1.00 










GLIMPSE II 


fixed [3.6] 


2.76 ±0.10 


1.66 ±0.03 


1.00 


0.49 ±0.01 


0.38 ±0.02 


0.28 ±0.03 


0.37 ±0.04 


fixed [4.5] 


2.75 ±0.15 


1.65 ±0.05 


1.00 


0.50 ±0.04 


0.39 ±0.01 


0.28 ±0.07 


0.37 ±0.08 


fixed [5.8] 


2.57 ±0.11 


1.58 ±0.03 


1.00 


0.55 ±0.01 


0.45 ±0.02 


0.36 ±0.01 


0.44 ±0.03 


fixed [8.0] 


2.61 ±0.11 


1.60 ±0.03 


1.00 


0.54 ±0.03 


0.43 ±0.03 


0.33 ±0.04 


0.42 ±0.01 


mean 


2.70±0.13 


1.63±0.03 


1.00 


0.53±0.02 


0.40±0.02 


0.32±0.04 


0.42±0.05 


Indebetouw 


2.50± 0.15 


1.55± 0.08 


1.00 


0.56± 0.06 


0.43 ± 0.08 


0.43±0.10 


0.43± 0.10 


Nishiyama 


2.86 ±0.08 


1.60 ±0.04 


1.00 


0.49±0.01 


0.39±0.01 


0.36±0.01 


0.42±0.01 


Gao 






1.00 


0.63± 0.01 


0.57± 0.03 


0.49± 0.03 


0.55± 0.03 



Notes. The center wavelength are those of the 2MASS filters. In bold face the assumed extinction A a /A Ks is indicated. 




10 5 -5 -to-0.2-0.1 o.o 0.1 0.2 



0.35 0.41 0.47 0.53 0.59 0.65 

Fig. 24. Galactic longitude vs. galactic latitude map of the A^j/A^ extinction coefficient assuming Aj/A Ks = 2.86 as long as longitudinal and 
latitudinal distribution of A A /A Ks (A equals to 3.6, 4.5, 5.8 and 8.0 yum), (a): The distribution of A A /A Ks as function of the galactic longitude. 
Different colors denotes different wavelengths, the zero line stands for fixed AjIA Ks and the black line indicates the data from Gao et al. (2009). 
(b): The map of the Ap 6i /A Ks extinction coefficient, x-axis denote galactic longitude while y-axis galactic latitude, (c): The same as (a) but for the 
galactic latitude. 



7. Conclusion 

Using an improved version of the Besancon model, we present 
here 3D extinction maps in the J, H, Ks, [3.6], [4.5], [5.8] 
and [5.8] bands using GLIMPSE II and VVV data. All 
the extinction maps are available online at CDS (Table 1 
and Table 2), as well as through the BEAMER webpage 
(http://mill.astro.puc.cl/BEAM/calculator.php). We derived new 
temperature -color relation for M giants which match better the 
observed color-magnitude diagrams. 

We present the 3D maps from 1 .25 //m until 8 fim. Due to 
the high sensitivity of the VVV data, we are able to trace large 
extinction until lOkpc. Our maps integrated along the line of 
sight up to 8 kpc show an excellent agreement to the 2D-maps 
from Schultheis et al. (2009) and Gonzalez et al. (2012). These 
maps show the same dust features and they show consistent Ak s 
values. 

Using the initial value from Nishiyama et al. (2009), we de- 
rived the mean extinction coefficient of our field in the seven 
bands, Aj/A Ks = 2.70±0.13, A H /A K , = 1.63+ 0.03, A [36] /A K , 
= 0.53+0.02, A l4 .5)/A Ks = 0.40±0.02, A [5 ,^/A Ks = 0.32+0.04, 
A[8d]/A fo = 0.42±0.05. The variation of the coefficients indi- 
cates that the wavelength dependence of interstellar extinction in 
the mid-IR varies from different lines of sight, which means that 



there is no "universal" IR extinction law. This is also in agree- 
ment with the study of Gao et al. (2009). 
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Appendix A: 3D reddening maps 
Appendix B: 2D reddening maps 
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Fig. A.l. 3d reddening and error maps for J-Ks VVV. X-axis denote galactic longitude while y-axis galactic latitude. The units are in <5E(J - 
Ks)kpc-'. 
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Fig. A.2. 3d reddening and error maps for E(H-Ks) VVV. 
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Fig. A.4. 3d reddening and error maps for E(Ks-[4.5]), 
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Fig. A.5. 3d reddening and error maps for E(Ks-[5.8]). 
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Fig. A.6. 3d reddening and error maps for E(Ks-[8.0]). 
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Fig. B.l. 2d reddening and error map for E(Ks-[3.6]), E(Ks-[4.5]), E(Ks-[5.8], E(K-[8.0]) integrated to 8 kpc distance. 
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